arXiv:1504.02107vl [astro-ph.CO] 8 Apr 2015 


Astronomy & Astrophysics manuscript no. Zevol © ESO 2015 

April 10, 2015 


The evolution of the spatially-resolved metal abundance 
in galaxy clusters up to z = 1.4 
(Research Note) 

S. Ettori 1 - 2 , A. Baldi 3 ,1. Balestra 4 , F. Gastaldello 5 ' 6 , S. Molendi 5 , P. Tozzi 7 

1 INAF, Osservatorio Astronomico di Bologna, via Ranzani 1,1-40127, Bologna, Italy 

2 INFN, Sezione di Bologna, viale Berti Pichat 6/2,1-40127 Bologna, Italy 

3 Physics and Astronomy Dept., Michigan State University, East Lansing, MI, 48824 USA 

4 INAF, Osservatorio Astronomico di Trieste, via G.B. Tiepolo 11,1-34131, Trieste, Italy 

5 INAF-IASF, via Bassini 15, 20133 Milan, Italy 

6 Department of Physics and Astronomy, University of California at Irvine. 4129 Frederick Reines Hall, Irvine, CA 92697^1575, 
USA 

7 INAF, Osservatorio Astrofisico di Arcetri, Largo Enrico Fermi 5,1-50125 Firenze, Italy 
Accepted 1 April 2015 


ABSTRACT 

Context. We present the combined analysis of the metal content of 83 objects in the redshift range 0.09-1.39, and spatially-resolved 
in the 3 bins (0-0.15,0.15-0.4, > 0.4)R 500 , as obtained with similar analysis using XMM-Newton data in Leccardi & Molendi (2008) 
and Baldi et al. (2012). 

Aims. By combining these two large datasets, we investigate the relations between abundance, temperature, radial position and redshift 
holding in the Intra-Cluster Medium. 

Methods. We fit functional forms to the combination of the different physical quantities of interest, i.e. ICM metal abundance, radius, 
and redshift. We use the pseudo-entropy ratio to separate the Cool-Core (CC) cluster population, where the central gas density tends 
to be relatively higher, cooler and more metal rich, from the Non-Cool-Core systems. 

Results. The average, redshift-independent, metal abundance measured in the 3 radial bins decrease moving outwards, with a mean 
metallicity in the core that is even 3 (two) times higher than the value of 0.16 times the solar abundance in Anders & Grevesse 
(1989) estimated at r > 0.4R 50 o in CC (NCC) objects. We find that the values of the emission-weighted metallicity are well-fitted 
by the relation Z(z) = Zo (1 + z)~ y at given radius. A significant scatter, intrinsic to the observed distribution and of the order of 
0.05 - 0.15, is observed below 0.4R 500 . The nominal best-fit value of y is significantly different from zero (> 3cr) in the inner cluster 
regions (y = 1.6 ± 0.2) and in CC clusters only. These results are confirmed also with a bootstrap analysis, which provides a still 
significant negative evolution in the core of CC systems (P > 99.9 per cent, when counting the number of random repetitions which 
provides y > 0). No redshift-evolution is observed when regions above the core (r > 0.15/?soo) are considered. A reasonable good 
fit of both the radial and redshift dependence is provided from the functional form Z(r,z) = Z 0 (1 + (r/0.15R5oo) 2 ) - ^ (1 + z)~ y , with 
(Zo,/3,y) = (0.83 ± 0.13,0.55 ± 0.07,1.7 ± 0.6) in CC clusters and (0.39 ± 0.04,0.37 ± 0.15,0.5 ± 0.5) for NCC systems. 
Conclusions. Our results represent the most extensive study of the spatially-resolved metal distribution in the cluster plasma as 
function of redshift. It defines the limits that numerical and analytic models describing the metal enrichment in the ICM have to meet. 
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1. Introduction 

The hot, thin X-ray emitting plasma in galaxy clusters (i.e. the 
intra-cluster medium ICM) is enriched with metals ejected form 
supernovae (SNe) explosions through subsequent episodes of 
star formation and subsequent diffusion through several mech¬ 
anisms, like, e.g., ram-pressure stripping, galactic winds, out¬ 
flows from Active Galactic Nuclei, galaxy-galaxy interactions 
(e.g. [Schindler & Diaferio|2008| . 

X-ray observations provide direct measurements of the metal 
abundance in the ICM, as well as their radial distribution and 
variation as a function of time. These measures represent the 
‘footprint’ of cosmic star formation history and are crucial to 
trace the effect of SN feedback on the ICM at different cosmic 
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epochs (e.g. |Ettori 2005 Calura et al.||2007 Cora et aT| [2008 
Fabjan et al.|2010f! 

Several studies have addressed the radial distributions of 
metals in the ICM at low redshift (e.g. |Finoguenov et al.|2000 


De Grandi & Molendi 200 1| Irwin & Bregman 2001; Tamura 


et al.| |2004[ |Baldi et al.||2007| |Leccardi & Molendi| |2008; 
Snowden et al. 2008). Few others have constrained the evolu¬ 
tion of the metal abundance at z <; 0.3 (e.g. Balestra et al. 2007; 
Maughan et al. 2008; Anderson et al. 2009; a statistical anal¬ 
ysis of the combination of these different datasets is presented 
in Andreon 2012) In Baldi et al. (2012), we have presented the 
XMM-Newton analysis of 39 galaxy clusters at 0.4 < z < 1.4, 
covering a temperature range of 2 % kT ^ 12.8 keV. We were 
able to resolve their abundance in 3 radial bins. 

In this work, we combine this dataset with the one presented 
Leccardi & Molendi (2008[), that includes 44 objects at z < 


0.31, with gas temperatures between 2.9 keV and 11.3 keV. The 
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analysis performed on the XMM-Newton data in the two samples 
is identical and can be statistically combined to probe the ICM 
abundance as function of radial position and redshift. 

We adopt a cosmological model with // (l = 70 km/s/Mpc, 
Q m = 0.3, and fl,\ = 0.7 throughout. Solar abundance values 
published in Anders & Grevessej ( 1989) > are adopted for refer¬ 
ence. Confidence intervals are quoted at lcr unless otherwise 
stated. 


2. Data analysis 

Details on the X-ray data reduction and analysis are presented 
in |Leccardi & Molendi| ([2008) and BaldietaL ( |2012| >, where the 
interested reader can also find a discussion on the systematic ef¬ 
fects that could affect our estimates of the ICM metallicity. We 
remind that most of these measurements, being for high temper¬ 
ature (kT > 2 keV) systems at medium/high redshifts (z > 0.1), 
are mainly based on the strength of the iron K line (at the rest- 
frame energies of 6.6-7.0 keV). The full dataset analysed here is 
available at http://pico.bo.astro.it/ ~settori/ 
abun/xmm. dat. We briefly mention here that a simultaneous 
spectral fit was performed on spectra accumulated from the two 
MOS detectors over the energy range 0.7 - 8.0 keV and using 
a Cash statistics applied to the source plus background counts, 
with a modelled background emission. The free parameters of 
the thermal model apec in XSPEC ( Arnaud|1996| > are tempera¬ 
ture, abundance, and normalization. Local absorption is fixed to 
the Galactic neutral hydrogen column density, as obtained from 
radio data ( jKalberla et al.||2005[ >, and the redshift to the value 
measured from optical spectroscopy. 

Because of the better statistics in lower redshift clusters. 


the X-ray spectral analysis of most objects in Leccardi & 
Molendi ( 2008) 1 was performed using a finer radial binning 
with respect to Baldi et al. ( 2012)>. Thus, w e had to degrade 
the spatial resolution of Leccardi & Molendi (|2008)l abundance 


profiles in order to match its resolution with Baldi et al. 
( 2012| > and obtain uniform abundance profiles in three radial 
bins: O-O.I 5 R 500 , 0.1 5-0.4/?5qo, and > O. 4 R 500 , where R 500 is 
estimated by computing a global gas temperature which does 
not include the core emission (< 0.15/Loo) and using it in a 
scaling relation (e.g. Vikhlinin et al. 2006; see Baldi et al. 2012 
for details). The spatial degradation was carried out first by 
interpolating the abundance profile (and its relative error) over 
the radial points of the observed surface brightness profile Sir) 
of each cluster. Then, the surface brightness profile was used 
as a weighting factor for the interpolated abundance profile, 
together with the errors on the abundance profile itself, so that, 
in each resulting bin, regions with higher surface brightness 
and smaller uncertainties on the measure of the abundance 
would have a larger weight in the computation of the average 
abundance. This scheme reproduces an emission-weighted 
profile of the metal abundance. Hydrodynamic al simulations 
have shown that an emission-weighted metallicity matches the 
original value in simulations at better than 5% (e.g. Rasia et al. 


( 2008| >). Following this approach, the average abundance <Z>, 
and the relative error cr <z> , in the R mm - R max bin was computed 


as: 
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where Z, and cr Z| . are the abundance profile and the error on the 
abundance profile, respectively, w, = 5,/cr| is the weighting 
factor, and the sum is done over the n radial bins falling in the 
radial range /? ml „ - R max in the abundance profile of each cluster 
as measured in|Leccardi & Molendi|(|20()8]). 


3. Results on the spatially-resolved abundance 
evolution 

This work present our final effort to use archived data to con¬ 
strain the evolution of the spatially-resolved ICM abundance. 
Our analysis separates between the central regions, where cool 
cores (if any) appear, and the rest of the cluster atmosphere. We 
require that the spectrum from which the metallicity is measured 
contains, at least, 300 net counts. A total of 221 data-points 
corresponding to the spatially-resolved spectroscopic measure¬ 
ments of the metal abundance distribution in 83 galaxy clusters 
in the redshift range 0.092-1.393 is considered. By estimating 
the pseudo-entropy ratio cr = ( Tq/T i) x (EMq/EM \)~ 1 ^, where 
Tq and T i are the temperatures measured in the r < 0.15/?soo 
region and in the 0.15-0.4i?5oo annulus, respectively, and EMq 
and EM\ are the corresponding emission measures, we define 
as non-cool-core (NCC) clusters the 54 (out of 83) objects that 
have cr > 0.6 (see e.g. Leccardi et al. (2010 1 and Baldi et ST) 
(2012b)). The data points are plotted as function as radius and 
redshift, separating between CC and NCC systems, in Fig. |T| 
All the results are presented in Table [T] 

In the r < 0.15/Coo radial bin, we have 70 data points over 
the entire redshift range. The distribution of the metal abun¬ 
dance correlates significantly with the redshift (Spearman’s rank 
p- = -0.25, corresponding to a significance of 2.lcr versus the 
null-hypothesis of no correlation; see Table |2j, with the tem¬ 
perature ( p T = -0.3, 2.5cr) and with the pseudo-entropy ratio 
(pc- = -0.4, 3.3cr). These correlations point towards a signif¬ 
icant effect induced from the cooling activities taking place in 
these regions, where the ICM at high density radiates more effi¬ 
ciently, lowering its global temperature and increasing its metal 
budget (see e.g. |Leccardi & Molendi| ( |2008] >). Fitting the data 
points in this bin with a power-law in the form Z oc (1 + z)~ 7 , we 
obtain y = 1.60 ± 0.22, different from zero at high significance 
(> 6cr). When a bootstrap analysis is done, i.e. the fit is repeated 
10 5 times after a random sampling with replacement and the me¬ 
dian, the 1st and 3rd qualifies are considered to estimate location 
and dispersion of the distribution of the best-fitting parameters, 
the central values are recovered but the relative errors are in¬ 
creased by about a factor of 2, mostly because a more proper 
sampling of the intrinsic scatter is performed with the bootstrap 
analysis. The negative evolution (y > 0) is still detected in the 
>99.9 per cent of the replica. When we consider CC and NCC 
systems separately, the normalisation Zq is different by a factor 
of 2.2 in favour of CC systems, where a negative evolution is 
still detected at very high significance. Instead, no evolution is 
measured (y = 0.5 ± 0.4) in NCC objects, although the case for 
negative evolution is still observed in about 80% of the boot¬ 
strapped repetitions. 

Eighty-three data points are located in the radial bin 
0.15/?5 oo < r < O. 4 R 500 , where a 2.8cr significant correla¬ 
tion between pseudo-entropy ratio and redshift is detected, sug¬ 
gesting an increase of cr (i.e larger incidence of NCC systems) 
with the redshift. A mild anti-correlation between Z and redshift 
is also present in this bin (significance of 2.2cr), We measure 
y = 0.70 ± 0.32 (lcr range of 0.16 - 1.31 after bootstrap anal¬ 
ysis). The case for negative evolution is still observed in about 
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Fig. 1 . Abundance distribution as function of redshift as estimated in 3 radial bins (O-O.I 5 R 500 , 0.15-0.4/^no, and > O. 4 R 500 ). Filled 
points indicate CC clusters. (Left) All data are plotted with the best-fits done in the given radial bin (dotted line) and including the 
radial dependence (dashed line); (right) the best-fits in redshift and radius (dashed lines) are shown with representative data points 
obtained as a weighted mean, with the relative error (solid bars) and dispersion (dashed bars). 


Table 1 . Results of the combined analysis on the redshift and radial distribution of the metal abundance. The columns show: the 
selected sample; the number of fitted data points; the total x 2 value; the intrinsic scatter estimated by requiring that the reduced x 2 
is equal to 1; the best-fit parameters of the adopted functional form Z(r,z ) = Zy (1 + (r/0.15/?5oo) 2 ) - ^ (1 + z)~ r . Median, 1st and 3rd 
quartiles from the distribution of the best-fitting parameters in the bootstrap analysis over a sample of 10 5 repetitions are indicated 
in the round parentheses. In the last column, labelled “constant”, we quote the weighted-mean, its relative error and the dispersion 
around it in the round parentheses. 


Sample 

N 

x L 

CT 

Zo 

7 


constant 

0 - 0.1 5 R 500 

70 

267.1 

0.151 ±0.020 

0.648 ±0.031(0.647 + o °“) 

1.60 ±0.22(1.59+°^) 

- 

0.451 ± 0.009(0.213) 

only CC 
only NCC 

27 

43 

140.7 

50.6 

0.159 ±0.031 
0.053 ± 0.038 

0.795 ± 0.046(0.780+°™) 
0.360 ±0.041(0.365!° ™) 

2.19 ± 0.28(2.10+°+i) 
0.45 ±0.40(0.51!°’!) 

- 

0.493 ±0.011(0.207) 
0.320 ±0.019(0.187) 

0.15 - 0.4R 5 oo 

83 

147.7 

0.044 ±0.012 

0.261 ± 0.022(0.271+°°!+) 

0.70 ± 0.32(0.66+°“) 

- 

0.220 ± 0.009(0.134) 

only CC 

29 

27.9 

0.012 ± 0.025 

0.287 ± 0.032(0.293!° °g) 

0.30 ± 0.47(0.33+1;!) 
0.37 ± 0.46(0.37!°“) 

- 

0.269 ±0.013(0.120) 

only NCC 

54 

96.8 

0.048 ± 0.017 

0.197 ±0.027(0.218!°;™) 

- 

0.178 ±0.012(0.146) 

> 0.4R 500 

68 

74.4 

0.029 ± 0.028 

0.168 ±0.028(0.173+° °!?) 

0.26 ±0.61(0.33+°!?) 
-0.88 ± 0.84(-0.71+ [■;!!) 

- 

0.158 ±0.014(0.170) 

only CC 

27 

40.1 

0.071 ±0.038 

0.131 ±0.031(0.138+™) 

- 

0.158 ±0.019(0.183) 

only NCC 

41 

30.7 

0.000 ± 0.008 

0.226 ± 0.053(0.223!°;™) 

1.39 ±0.92(1.31!°;“) 

- 

0.157 ±0.020(0.164) 

all 

221 

506.3 

0.092 ± 0.010 

0.702 ±0.031(0.699+° °!') 
0.838 ± 0.046(0.831+°!!,) 

1.33 ± 0.18(1.31+° !i) 
1.75 ± 0.24(1.72+°") 

0.56 ±0.03(0.55+°!!!) 

0.300 ±0.006(0.181) 

only CC 

83 

225.1 

0.105 ±0.015 

0.55 ± 0.04(0.55+°!!!) 

0.365 ± 0.008(0.200) 

only NCC 

138 

185.3 

0.037 ±0.013 

0.387 ± 0.038(0.393+°™) 

0.52 ± 0.29(0.53+°!+) 

0.39 ± 0.06(0.37+°!+) 

0.207 ±0.009(0.174) 


91% of the bootstrapped samples. Similar values are measured 
in the population of NCC clusters, which dominates the sample 
in this radial bin, with a normalisation (and a mean metallicity 
value) that is lower by about 30 per cent than the corresponding 
value measured for CC systems. 


Finally, we do not find evidence of negative redshift evolu¬ 
tion of the metallicity also in the last radial bin considered in 
our analysis (r > 0.4/Coo), where we constrain a value for y of 
0.26 ± 0.61 (y > 0 with a probability P ~ 65% after bootstrap 
analysis) with 68 data points. Similar results are obtained in the 


CC and NCC populations, where an average value of Z = 0.16Z 0 
is estimated. 


Being the sample highly heterogeneous and not representa¬ 
tive in any sense of the cluster population, it is not expected to 
weight properly for the relative incidence of these cool cores 
clusters. Moving at higher redshift, for instance, lower metal 
abundance associated to the inner regions suggest that less pro¬ 


nounced cool cores are resolved in our sample (see e.g. Santos 
[eTaLl ( |20T0l i; |Samuele et aL] ( |20TT] i) 
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Table 2. The Spearman’s rank correlation coefficient (and 
the significance in terms of the standard deviations by which 
the sum-squared difference of ranks deviates from its null- 
hypothesis expected value) between metallicity and redshift (p z ), 
metallicity and temperature (pr), metallicity and pseudo-entropy 
ratio Ov). and pseudo-entropy ratio and redshift (p, r ,), respec¬ 
tively. The pseudo-entropy ratio is larger in NCC systems. 


Sample p z p T p^ p lnz 

0-0.15 R 500 -0.25 (2.07) -0.31 (2.54) -0.39 (3.27) 0.18 (-1.51) 

0.15 - 0.4/? 50 o -0.24(2.18) -0.13 (1.14) -0.13 (1.18) 0.31 (-2.80) 

> 0.4f?5o 0 -0.20 (1.63) -0.09 (0.75) -0.10 (0.84) 0.23 (-1.85) 

all -0.20(3.04) -0.09 (1.27) -0.20(2.93) 0.25 (-3.66) 


4. Summary and Conclusions 

We present the statistical analysis of the spatially resolved spec¬ 
troscopic results of the XMM-Newton observations of 83 clusters 
of galaxies in the redshift range 0.09 < z < 1-39 and with a gas 
temperature between 2 keV and 12.8 keV. In a consistent and 
similar way (see details in |Leccardi & M olendi (2008) and Baldi 
|et al.| ( |20T2| >), the cluster emission has been resolved for all the 
objects in the dataset in three different regions: the core region 
(corresponding to 0 < r < 0.1 5/Goo). the region immediately 
surrounding the core (0.15/Goo < r < 0.4/GooX and the outskirts 
of the cluster (r > 0.4/Goo)- Using the pseudo-entropy ratio cr, 
we also distinguish between Cool-Core and Non-Cool-Core ob¬ 
jects, with CC clusters that represent -35 per cent of the whole 
sample. Our results are summarised in Tables [I] and [2] 

We find strong (significance at 3<x or higher) correlations be¬ 
tween (i) metallicity and pseudo-entropy ratio (in the inner bin) 
and, in the entire sample, (ii) metallicity and redshift, and (iii) 
pseudo-entropy ratio and redshift. 

By fitting a constant value (under the assumption of no¬ 
evolution) to the metal abundance measured in the 3 radial bins, 
we confirm a neat decrease moving outwards, with the inner bin 
being even 3 (two) times more metal-rich than the outermost bin 
in CC (NCC) systems (0.49/0.32, 0.27/0.18 and 0.16/0.16 for 
CC/NCC systems in the 3 radial bins, respectively). 

We constrain the parameters of the functional form Z(r,z ) = 
Zo (1 + (r/0.15f?5oo) 2 )”^ (1 + z)~ y . The best-fitting values are 
Zo = 0.70 ± 0.03 (0.62 - 0.79, lst-3rd quartile from the boot¬ 
strap analysis), /? = 0.56 ± 0.03 (0.48 - 0.65) and y = 1.3 ± 0.2 
(0.95 - 1.72). CC and NCC systems behave differently, with 
(Zo,y8,y) equals to (0.83 + 0.13,0.55 + 0.07,1.7 + 0.6) and 
(0.39 ± 0.04,0.37 + 0.15,0.5 + 0.5), respectively (error bars refer 
here to the bootstrap analysis). 

At a given radius, we find a statistically significant negative 
evolution in the distribution of the metal abundance as function 
of redshift in the inner radial bin (> 3.5cr; the probability from 
the bootstrap analysis to measure y > 0 is larger than 99.9 per 
cent). We show that this level of evolution with the cosmic time 
is entirely consistent with the same evolution observed in the 
CC-only population. On average, we measure a mean metallicity 
in the local systems that is about 3—4.5 (1.4 in NCC clusters) 
times the value measured at z ~ 1, which is consistent with the 


constraints we publish in Baldi et al. 

2012) and, using (mostly) 

Chandra data, in Balestra et al. (2007 

). At the same redshift, the 


metallicity decreases radially, reaching at ~ 0.5Rson about 37 (in 
CC) and 50 (in NCC) per cent of the value estimated at 0.15/Goo- 
By requiring a reduced^ 2 of 1 in the performed fit, we estimate 
an intrinsic scatter in Z of about 0.15,0.05 and 0.05 in the 3 
radial bins, respectively. In the inner bin, the scatter of ~ 0.15 


is measured also in the CC systems only, whereas NCC clusters 
show a value of 0.05 consistently in the inner 2 radial bins. 

CC clusters show evidence of negative evolution just in the 
inner radial bin (< 0.15/Goo). where the accumulation of the 
plasma during the cluster formation makes the ICM denser and, 
as consequence of radiative processes, cooler and more metal- 
rich. Most of the above-mentioned correlations (i.e. p z , p-y and 
Per in Tab. [2]i can be then interpreted as due to the strong in¬ 
fluence of this centrally-located cool gas. NCC systems do not 
show any statistically significant deviation from y - 0 at any ra¬ 
dial bin. Only in the outer bin (r > 0.4/Goo). both CC and NCC 
systems present similar average metallicity, and an intrinsic scat¬ 
ter almost consistent with zero. We conclude that, for what con¬ 
cerns the metal distribution, this cluster volume is the most ho¬ 
mogeneous in both CC and NCC systems, and might be consid¬ 
ered the most representative of the average cluster environment. 

We stress that the constraints presented in this work repre¬ 
sent the most robust ones reachable with archived XMM-Newton 
and Chandra data on the distribution of the metal abundance as 
function of radius and redshift in the ICM. Larger samples of 
high redshift X-ray clusters, together with deeper Chandra and 
XMM-Newton observations even of the clusters already present 
in the archives, would be crucial to put tighter constraints on 
the evolution of the ICM metal abundance and to provide a ro¬ 
bust modelling of the physical processes responsible for the en¬ 
richment of the cluster plasma during its assembly over cosmic 
time, from the accumulation of the iron mass in the CC cluster 
cores (e.g. De Grandi et al. 2004, Bohringer et al. 2004) to the 
history of the processes responsible for the metal release into 
the ICM through ejections from supernovae type la and II (e.g. 
Ettori 2005) and through galaxy transformation in the cluster 
volume (e.g. Calura et al. 2007), and their modelling in joint 
semi-analytic/hydrodynamical numerical simulations (e.g. Cora 
et al. 2008, Fabjan et al. 2010). 
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